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Abstract 

> 

("^ , An order parameter description of the Anderson-Mott transition (AMT) is 

I given. We first derive an order parameter field theory for the AMT, and 

■r^-j- ■ then present a mean-field solution. It is shown that the mean-field critical 

0^ ■ exponents are exact above the upper critical dimension. Renormalization 

group methods are then used to show that a random-field like term is gen- 
erated under renormalization. This leads to similarities between the AMT 
^ ^ _ and random- field magnets, and to an upper critical dimension = 6 for the 

■ AMT. For d<6, ane = 6 — d expansion is used to calculate the critical 

O ■ exponents. To first order in e they are found to coincide with the exponents 

for the random-field Ising model. We then discuss a general scaling theory for 
the AMT. Some well established scaling relations, such as Wegner's scaling 
^ ' law, are found to be modified due to random- field effects. New experiments 

' are proposed to test for random-field aspects of the AMT. 

PACS numbers: 71.30.-hh, 64.60.Ak, 75.10.Nr 
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I. INTRODUCTION 



Metal-insulator transitions of purely electronic origin, i.e. those for which the structure 
of the ionic background does not play an essential role, are commonly divided into two 
categories. In one category the transition is triggered by electronic correlations, in the other 
it is driven by disorder.^ The first case is known as a Mott transition, and the second one 
as an Anderson transition. It is believed that for many real metal-insulator transitions both 
correlations and disorder are relevant, and this expectation is also borne out by model studies 
of disordered interacting electron systems.il The resulting quantum phase transition, which 
carries aspects of both types of transitions, we call an Anderson- Mott transition (AMT). 

Most of the theoretical work to date on the AMT has been based on a small disorder 
expansion near d = 2, which is believed to be the lower critical dimension for the AMT.i A 
generalized matrix nonlinear a-model (NLaM) has proven very useful for this purposeJl This 
model is a generalization of Wegner's theory for the Anderson transition in noninteracting 
electronic systems.^ It has led to the identification of various universality classes for the 
AMT, as well as to an e = d — 2 expansion for the critical exponents.!!! Complementary 
work has recently been done by adding disorder to models which display a Mott transition. 
In the limit of infinite dimensions the resulting models have been studied in some detail.i 

From a standard phase transition point of view these approaches to the AMT are quite 
unconventional. The standard procedure for solving phase transition problems! is to (1) 
identify the relevant order parameter (OP), and possibly derive or postulate an effective 
field theory for the OP, (2) construct a Landau, or mean-field description of the phase 
transition, and (3) use renormalization group (RG) methods to identify the upper critical 
dimension, ci^, and to compute the critical exponents in an e = — d expansion. RG 
methods can also be used to derive a general scaling description of the phase transition for 
d<d+. 

Early attempts to implement this standard program for the Anderson transition! failed, 
because the most obvious simple OP, viz. the single-particle density of states (DOS) at the 
Fermi level,0 turned out to be uncritical.0 As a result, the Anderson transition can onlv 
be described by the NLcrM approach, or possibly in terms of a complicated functional 0P.liHl 
The situation is, however, fundamentally different at the AMT, where the DOS is generally 
believed to be critical, and to vanish at the transition. This raises the prospect of using the 
DOS for a conventional OP description of the AMT, making the AMT conceptually simpler 
than the Anderson transition. 

In two previous short publicationsEl'0 we have shown that these expectations are indeed 
justified, and have discussed such an OP approach to the AMT. The purpose of the present 
paper is to explain the technical details of the OP approach, and to present a number of 
additional results. 

One of the most far-reaching implications of our approach is that the AMT is in some 
respects similar to magnetic transitions in random fields.0 We will present the technical 
reasons underlying this conclusion in Sec. |V| below. Here we argue that a random-field 
structure of the theory could have been anticipated on general physical grounds. Let us 
consider a model of an interacting disordered electron gas. In terms of anticommuting 
Grassmann fields, ijj and ifj, the action can be written Jl 
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Here x = (x, r) with r denoting imaginary time, J dx = / dx Jq dr, m is the electron 
mass, fi is the chemical potential, cr is a spin label, and for simplicity we have assumed 
an instantaneous point-like electron-electron interaction with strength F. m(x) is a random 
potential which represents the disorder. We assume u to be (5-correlated, and to obey a 
Gaussian distribution with second moment 

HxMy)} = — ^5(x-y) , (1.1b) 

where the braces denote the disorder average, Np is the bare DOS per spin at the Fermi 
energy, and Tei is the bare elastic mean-free time. For future reference we denote the disorder 
related part of the action by 

Sdis = - Y, dx m(x) iIj^{x)4j^{x) = - ^ / rfx u(x) ?/'^,„(x)?/'^,„(x) . (1.2) 

In the second equality in Eq. ( |1.2|) a Matsubara frequency decomposition of V'(r) and '?/'(r) 
has been used. 

As mentioned above, the most obvious OP for the AMT is the single-particle DOS, 
A^, at the Fermi level. In terms of Grassmann variables this quantity is proportional to 
the zero- frequency limit of the expectation value of the composite fermionic variable ipip: 
N = lmN{iuJn + with 



where we have normalized the DOS by 2Np. Equations (|1.2| , |1.3|) suggest that the OP for 
the AMT couples directly to the random potential u, and that this random field (RF) term 
is structurally identical to the one that appears in magnetic RF problems. S'tl Notice that 
this structure appears in both interacting and noninteracting localization theories. However, 
in the interacting case there is an additional physical feature: The interaction will in general 
favor a local electron arrangement that is different from the one favored by the random 
potential. This kind of frustration is crucial for the generation of the typical RF effects 
which are known from RF magnets. One should therefore expect RF effects at the AMT, 
but not necessarily at the Anderson transition. 

This conclusion has several important implications. For instance, one expects hyper- 
scaling to be violated at the AMT due to a dangerous irrelevant variable, as it is in RF 
magnets.li^ We will see below that this is indeed the case and that, as a result, Wegner's 
scaling law0 relating the conductivity exponent, s, to the correlation length exponent, v, is 
modified. 

The plan of this paper is as follows. In Sec. || we review the NLcjM approach to the 
AMT, and use it to derive an OP field theory for the transition. In Sec. |IT1| the OP field 
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theory is solved in the saddle point or mean-field approximation. We show that the mean- 
field solution is stable for d > d'^, and that the mean- field exponents are exact in this regime. 
In Sec. |V| RG methods are used to prove that the theory has RF aspects, to establish that 
d'^ = 6, and to set up an e = d — 6 expansion for the critical exponents. In Sec. we discuss 
a general scaling description of the AMT. We conclude in Sec. [V| with a summary of the 
paper and a discussion of some open problems as well as of experimental implications of our 
results. 



II. FORMALISM 



In this section we start by briefiy reviewing the field theoretic description of disordered 
interacting electron systems. We focus on the NLcxM, which is a free field theory with 
nonlinear constraints relating massive and massless modes in the ordered phase, which in 
the case of the AMT is the metallic one. We then integrate out the massless excitations 
to obtain a field theory solely in terms of the massive modes, which are directly related to 
the OP for the transition, i.e. the DOS. This is different from the usual treatment of the 
NLo"M, which integrates out the massive fluctuations to obtain a theory for the soft modes. 
Our procedure is closely analogous to the treatment of the 0{n) symmetric NLcxM in the 
limit of large n.lli 



A. The Model 



We will start from the action given by Eq. (1.1a). The NLcrM for the AMT can be derived 



from Eq. ( |l.la]) by assuming that all of the relevant physics can be expressed in terms of 
long-wavelength and low-frequency fluctuations of the number density, the spin density, and 
the single-particle spectral density. Technically this is achieved by making long-wavelength 
approximations, and by introducing classical composite variables that are related to the 
operators mentioned above. The quenched disorder is handled by means of the replica trick. 
The resulting action reads,ili 

1 r . /_^, .\2 



S[Q] = -— J dx tr(vg(x))' + 2H Jd^ tr{QQ{^. 

jd^ [Q(x)7(")Q(x)] , (2.1a) 



ttT 



u=s,t 

where 

[Q(x)7(^)g(x)] =K, 5n.+„3,n.+n, E E l")'^ 

n\,n2,n-i,n4, a r=0,3 

^tr[{Tr®so)Ql%,)tr{{Tr®s,)QT,n,) , (2.1b) 

and 

[Q(X)7(*)Q(X)] = -K, E ^n,+n3,n.+n, E E E 

n-i,n2,n-j,,rn a r=0,3 i=l 



xtr((r,®s,)Q^^„Jtr((7:,®.,)Q^^„J , (2.1c) 
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Here Q is a classical field that is, roughly speaking, composed of two fermionic fields. It 
carries two Matsubara frequency labels, n and m, and two replica labels, a and 13. The matrix 
elements Q^m ^Pi^ quaternions, with the quaternion degrees of freedom describing the 
particle-hole {Q ~ iIjiIj) and particle-particle {Q ~ ifj^l)) channels, respectively. We will 
restrict ourselves to the particle-hole degrees of freedom. The matrix elements can then be 
expanded in a restricted spin-quaternion basis. 



3 



= E E IQli (rr ® s,) , (2.2) 



r=0,3 1=0 



with To,i,2,3 the quaternion basis, and So,i,2,3 the spin basis (si^2,3 are the Pauli matrices). 
The matrix Q is subject to the nonlinear constraints, 

Q2 = 1 , (2.3a) 
tr Q = , (2.3b) 

gt = c^qTc = Q . (2.3c) 
The last equation expresses the requirements of hermiticity and charge conjugation; the 



matrix C is block-diagonal with elements iri ® S2- Equation ( 2.3c ) implies. 



Qnt = i-r srrQ:i , (2.4b) 



with 5*0 = 1 and 5*1^2,3 = — 1- 

In Eqs. (|2.1| ), G = 2/'7ia, with a the bare conductivity, is a measure of the disorder, and 
H = 7tNf/2 is a frequency coupling parameter. Kg and Kt are bare interaction amplitudes in 
the spin singlet and spin triplet channels, respectively, and = Snm^af^ (tq ® sq) ujn, with 
ujn = 27rTn, is a bosonic frequency matrix. Notice that Kg < for repulsive interactions. 

The correlation functions of the Q determine the physical quantities. Correlations of Qnm 
with nm < determine the diffusive modes which describe charge, spin, and heat diffusion, 
while the DOS is determined by (Q""), cf. Eq. ( ^.3aD below. It is therefore convenient to 
separate Q into blocks, 

= e(nm)QKx) + e{n)Q{-m)qfM 

+e(-n)e(m)(gt)^^(x) . (2.5) 

Normally a NLaM is treated by integrating out the massive modes, i.e. the Qnm in this case, 
to obtain an effective theory for the massless modes, which here are the diffusion processes 
described by q and However, since our goal is to obtain a field theory for the OP for the 
AMT, Qnn, "we will proceed differently, and instead integrate out the massless g-fields. To 
this end, we use a functional integral representation of the delta function, 
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n m\^) - 1] = / ^[A] exp{-^ / rfx tr (A(x)[g2(x) - 1 



(2.6) 



where tr denotes a trace over all discrete indices, and the factor of 1/2G has been inserted 
for convenience. Together with Eq. (|2.1a|) this allows us to write the action as, 



5i[Q,A] = — /rfxtr 



A(x)[g2(x)-ii] + (9,g(x 



ttT 



+ 2H J dx tr[nQ{x] 

y: [q(x)7(")q(x)] 



(2.7) 



u=s,t 



We note that the constraint given by Eq. (|2.3b| ) does not involve q, so it can be imposed 
later.Ei Furthermore, the contraints expressed by Eqs. ( |2.3c| ) or ( [^.4D do not couple indepen- 
dent blocks of Q. Finally, by decomposing A into blocks like Q one sees that the elements 



Anm with nm > 0, together with the tracelessness condition, Eq. ( 2.3b ), are sufficient for 
enforcing the constraint Q"^ = 1. We can therefore restrict ourselves to Anm with nm > 0. 
The partition function of the replicated theory (N replicas) is given by. 



= I D[Q] exp{S[Q]) 
D[Q] D[A] exp(5i[g,A]) = / D[Q] D[A] exp(52[g,A]) . 



(2.8) 



(2.9) 



by integrating out the massless g-fields. Since the action 5*1, Eq. (|2.71 ), is quadratic in these 
fields, this can be done exactly. Formally one obtains. 



In the last equality in Eq. (p.8|) we have defined yet another action, 

S2[Q, A] = \n jD[q,q^ exp(^i[g,A]) , 



S2[Q,A] = S,[Q,A]-^Tr lnM[A] 



(2.10a) 



Here Tr denotes a trace over all degrees of freedom, and M is a matrix with elements, 
--M--fi- fx) = ^--^ Y: tr {rr,Tr,T}j tr (.,3.,4) -A--(x) 



rir2 nin2,n3n4 



2G 



2G 



L+n4,n2+n3 



2'kTGK^. 



;2.iob) 



with = s and 1^1,2,3 = t. In giving Eq. (|2.10a|) we have omitted some terms which result 
from the interaction, i.e. the last term in Eq. ( p^.7| ), coupling Q and q. These terms will turn 
out to be irrelevant (in the RG sense) because they are of higher order in the frequency. 



6 



B. The Order Parameter Field Theory 



The field theory given by Eqs. (|2.10|) has been formulated in terms of the OP for the 
AMT and the auxiliary, or ghost, field A. In order to obtain a theory solely in terms of Q we 
need to integrate out A. We will do this perturbatively, and verify later that the neglected 
terms are either RG irrelevant, or provide analytic corrections to nonuniversal coefficients 
in the theory. Either way they do not modify the critical behavior. Physically this is to be 
expected, since it is easily established that the A-modes are not critical near the AMT.0 

To proceed we write the average of A as 

CA°^(x)) = 5,o5.o5„„.C , (2.11a) 
with In to be determined later. We write 

;aKx) = (;aKx)) + ;C^(x) , (2.11b) 
and expand in powers of ip to determine a new OP action, 6*3 [Q], which is defined by 

^3[Q] = In 1 D[i;] exp(S2[g, (A) + ^j]) . (2.12) 

Integrating over ip in Gaussian approximation, neglecting constant contributions to the 
action, and neglecting explicit interaction terms (see the end of this subsection for a justifi- 
cation), we obtain 

Ss[Q] = J rfx tr [(VQ(x))' + (A)g2(x)] +2Hjd^ tr [nQ(x)] 

I tr [(1 - /)g2(x)] - ^ Jd^tr g^(x) 

Jd^it^+Q'i^))it^-Q'i^)) > (2-13) 



u 

4 



2G2 



where tr± denotes 'half-traces' that sum over all replica labels but only over positive and 
negative frequencies, respectively: tr+ = J2a J2n>o 5 t^- = J2a J2n<o ■ f = /((^)) is a 
matrix with elements If^mi^) = ^roko^nmfn with. 



4ip„4l,[p2 + i(C + C)]^ .=0 



ni=0 ~^ 2^^"! '^n-i-n+m) 



(2.14) 



for n > 0, where /p = (27r)^'^/ rfp. For n < 0, the same expression holds with the sum over 
m from to cxd, and with n and m interchanged in the summand. The coefficients u and v 
in Eq. ( |2.13| ) are given by, 

M = -G'/(d/„/d4)„=o ' (2-15a) 

and 

v = -u' f 1 . (2.15b) 
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In writing Eq. (|2.13|) we have localized the coefficients u and v in space and imaginary 
time. Frequency and wavenumber dependent corrections to this result turn out to either 
renormalize nonuniversal coefficients, or to be RG irrelevant. Note that v is infinite for 
d < 4. u also contains singularities in d < 4. In the main text of this paper we restrict 
ourselves to the region d ^ d'^ = 6, and ignore these singularities. In Appendix ^ we show 
that, at the mean- field level, they are of no consequence for the critical behavior as long as 
d > 3. Beyond the mean-field level, these singularities are one of several severe difficulties 
one encounters. This point will be further discussed in Sec. Finally, we note that our 
integrating out of A has obscured the condition 

W = , (2.16a) 



which follows from Eq. ( p.llb|) . By returning to the action 5*2 one sees that Eq. (p. 16a]) 
implies, 

= I - - S^ {IQfMYrQili^)) + . . . , (2.16b) 

I3,m,r,i 



with Si given after Eq. ( |2.4b| ). Corrections to Eq. ( p.l6b| ) come from higher orders in the 



■j/'-expansion. We will show below that they do not change the critical behavior. 

We conclude this section by mentioning two approximations that were used to derive 
Eqs. ( p.l3| ) and ( p.l6| ). First, the coefficients given in Eq. ( 2.13|) follow by expanding 5*2 in 



Eq. ( |2.12D to 0{ip'^), neglecting all terms of higher order in ip. Second, all explicit interaction 
terms, such as the last terms in Eq. ( p.la| ), or the terms mentioned after Eq. (|2.10b| ), have 
been neglected in Eq. ( p.l3| ). In the next section we will use RG methods to argue that 
these approximations are justified, at least in the vicinity of d^. 



III. MEAN-FIELD THEORY 

In this section we solve the field theory derived in Sec. || in the mean-field or saddle 
point approximation. We will show that the mean-field theory describes a phase transition, 
and that it is an AMT. We then show that the saddle point solution is locally stable, and 
that the mean-field exponents are exact for d > d^. 



A. Mean-Field Solution 



We can use either Eqs. ( p.lOj ) or ( p.l3| ) - (|2.16|) to construct a mean-field theory. We 



choose to use the latter. We look for saddle point (SP) solutions Qsp, Agp that are spatially 
uniform and satisfy 



(O )°'^ 



(3.1a) 



i.{K)ni = SroSioSnmSa(,i^^'> , (3.1b) 

where the superscript (0) denotes the SP approximation. The replica, frequency, and spin- 
quaternion structures in Eqs. (|3.1|) are motivated by the fact that (^Q^m) and (* A^^) have 
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these properties (cf. Eq. ( p.llap ), and that in the mean-field approximation averages are 
replaced by the corresponding SP values. 
Equations ( |3.1| ) and ( |2.16b| ) yield 

(iVl°^f = l-/n((A)) . (3.2a) 
Taking the extremum of the action 5*3, Eq. (|2.13D, with respect to Q or iY*^") gives 



4°) = 2G/7a;„/iV(0) , (3.2b) 



where we have used Eq. ( p. 2a]) . 



Let us discuss some aspects of these results. First, the DOS at a frequency or energy fi, 
measured from the Fermi surface, and normalized by 2Nf (cf. Eq. (|1.3| )), is given by,u 

iV(fi) = Re (°Q^^(x)),^„^o+io = Re N^V^^^n+io ■ (3.3a) 
In mean-field approximation this yields 

iV(J])=iVW|,,„^f,+,o • (3.3b) 

Clearly, the DOS decreases with increasing disorder G. Furthermore, N{VL) is nonanalytic 
at i7 = 0, see Appendix 0. If we insert Eqs. (|3.2| ) into Eq. ( p.3b|) , and iterate to first order 
in G, then we just recover the well-known 'Coulomb anomaly' of the DOS.0 With further 
increasing disorder, N'^'^\VL = 0) eventually vanishes at a critical value Gc of G. The mean- 
field approximation thus describes a phase transition with a vanishing DOS at the Fermi 
level, the hallmark of the AMT. (Transport properties we will discuss shortly). Note that in 
the absence of interactions the quantity /, Eq. (|2.14|) , vanishes, and hence N^^^ = 1. This 
refiects the fact that the DOS does not vanish at an Anderson transitionlHI 

Second, it is easy to determine the behavior of N^^^ in the vicinity of Gc- This is true 
even though Eqs. ( |3.2D and ( |2.14| ) represent an integral equation for The point is that 
£rn in Eq. ( p.l4|) gets integrated over all frequencies. At nonzero m it is given by Um times 
a finite, nonuniversal prefactor which depends on A^^-*, see Eq. (|3.2b|) . This implies that we 



can treat im as a finite, nonuniversal quantity. The distance from the critical point, t, is 
proportional to G — Gc- In mean-field approximation t is given by, 

t(°) = l-/„=o . (3.4) 

Near the critical point we have 

Ar£ = (t(°))^/^ , (3.5) 
which yields the mean-field value of the critical exponent (3, 

(3=1/2 . (3.6) 



Note that the sign of A^^=o Eq. (3^) is determined by the way the zero- frequency limit is 
taken: A^„ is an odd function of n. In this respect the external Matsubara frequency acts like 
a symmetry breaking external field at the mean-field AMT. The DOS, of course, is positive 
definite since in order to obtain it from the branch cut must always be approached from 
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above, Eqs. (p.3|). To put it another way, Nn is i times a causal function that is an odd 
function of complex frequency (viz. the Green function). N{VL) is the spectrum of that 
function, which is even in Vt. 

Third, it is not much more difficult to obtain all other critical exponents. Let us write 
Q as its expectation value plus fluctuations, as we did with A in Eq. (|2.11b| ). 



where the factor of V2G has been inserted for later convenience. The action governing 
Gaussian fluctuations about the mean-fleld solution in the critical region then follows from 
Eq. (|]T3D as. 



t/x tr [(V0(x))' + ^ _ i{0)^2(^) 



V 



/rfx(tr+0(x))(tr_0(x)) + O(0^) . (3.8) 



All remaining critical exponents can now be read off Eq. ( p.8|) : Comparing the flrst and the 
third terms on the r.h.s. yields the correlation length exponent u = 1/2. With i^^^ u/Q, 
the flrst and the second term give the dynamical exponent z = 3. (This result one can also 
conflrm by explicitly calculating N{Q ^ 0) at G = Gc from the mean-fleld equation, see 
Appendix Finally, inspection shows that the 0-0 correlation function near the transition 
in the limit of small frequencies and small wavenumbers has a standard Ornstein-Zernike 
form, which yields exponents 7 = 1 and rj = 0. We thus have standard mean-fleld values for 
all static exponents, 

(3 = u=l/2 , 7=1 , 77 = , 5 = 3 , (3.9a) 

and for the dynamical exponent we have, 

z = 3 . (3.9b) 

Inspection of Eq. ( p.8|) further shows that the AMT saddle point is a local minimum and 
therefore stable. 

Apart from the quantities whose critical properties are governed by the exponents given 
in Eqs. ( |3.9|) we are also interested in the transport properties, and in some other ther- 
modynamic quantities. Let us flrst consider the charge diffusion coefficient, Dc. It can be 
obtained from the NLcrM by a direct calculation of the particle-hole spin-singlet q-q corre- 
lation function.! To this end we consider the action ^i, Eq. (|2.7|). If we replace Q and A by 
their mean-fleld values, Eqs. ( |3.1|) , ( p.2| ), we get an action that is quadratic in q. Inversion 
of the corresponding matrix yields the q-q propagator in mean-fleld approximation. It has 
a diffusive structure with diffusion constant 



r(0) 

Dc= .^r;.,(o) ' (3.10a) 



GH + GKsKlo 



that is, Dc vanishes like the OP. The same argument applied to the particle-hole spin-triplet 
q-q propagator and to a q-q correlation that is off-diagonal in replica space, which contain the 
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spin and heat diffusion coefficients Ds and Dh, respectively,! establish that these transport 
coefficients also vanish like the OP, 

GH + GKtNi% 

Dj, = Nr,=o/GH . (3.10c) 

Finally, we note that the SP approximation is not in any way related to a systematic 
small disorder (G) expansion. As a consequence, it is in general not possible to relate a G- 
expansion of the SP theory to previous work on the NLcxM^i, or to many-body perturbation 
theory.! An exception is the G-expansion of the OP itself, see the discussion after Eq. (|3.3b|) , 
and Eq. (^). 

B. The Mean-Field Transition as a Renormalization Group Fixed Point 

Here we discuss the mean-field solution presented in the last subsection in terms of a RG 
fixed point (FP). We will see that the mean- field critical behavior is exact for d > . We 
will also use the RG to derive scaling arguments which allow us to determine the critical 
behavior of thermodynamic susceptibilities, such as the specific heat coefficient and the 
density susceptibility dn/d^, and of the electrical conductivity. 

Consider the action given by Eq. ( |2.13|) , or Eq. ( |3.8D . Our parameter space is spanned 



by = {c,t,H,Ks^t,ui,U2}- Here c, whose bare value is 1/2G, is the coefficient of the 
gradient squared term in Eq. (|2.13| ), and t is a measure of the distance from criticality. Its 
bare value is t'^^^u/2G'^ with t^^^ given by Eq. (p^. ui = u/AG'^ and U2 = v/AG'^ are the 
quartic coupling constants, H is the frequency coupling, and Kg t are the interaction coupling 
constants. Since, in anticipation of the RG arguments below, we have dropped the explicit 
interaction terms from Eq. ( p.l3| ), this action contains Kg t only implicitly via the matrix 
f. The explicit interaction dependence is easily restored by adding to Eq. ( ^.13| ) the last 
term in Eq. ( |2.1a|) with Q replaced by Q. The contributions omitted from Eqs. ( p.lO|) yield 
a term of the same structure. Structurally, Eq. (|2.13|) is just a matrix version of standard 
0^-theory. Consequently, in looking for a mean-field FP, we follow Ref. M and |1^. We define 



the scale dimension of a length to be —1, and require the exponent 77, i.e. the anomalous 
dimension of the Q-field, to be zero. The exact scale dimension of Q is then [Q] = {d — 2)/2. 
Power counting then yields the scale dimensions of the various coupling constants, 

[c] = , [t] = 2 , [u,] = [u2]=A-d , 

[HT] = {d + 2)/2 , [Ks,tT] = 2 . (3.11a) 

Sometimes the dimension of temperature, T, or frequency, Q, is chosen to be [T] = [Q] = d, 
which results in [H] = 1 — d/2 and [Kg^t] = 2 — d. This choice lumps the anomalous part 
of the dynamical scaling exponent into the scaling behavior of H, and is convenient in 
considerations near the lower critical dimension (i^ = 2.i Near the upper critical dimension, 
d'^, the scaling dimension of H or its equivalent is usually chosen to be zero. Here we adopt 
this choice. Then 
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[T] = [n]=z = d/2 + l 



(3.11b) 



Here 5 is a 'naive' dynamical scaling exponent. As we will see below, the effective dynamical 
exponent z given in Eq. ( |3.9b|) results from z if the existence of dangerous irrelevant variables 
is taken into account. Equations ( |3.11aD ,( |3.11b| ) give 

= (2 - rf)/2 . (3.11c) 

With either choice for [T], [Ks^t] < for d > 2, that is, the electron-electron interaction 
is irrelevant. This justifies, a posteriori, our omitting the interaction terms in Eq. ( p. 131) . 



The same reasoning implies that the quartic coupling constants Ui^2 are irrelevant for d > A. 
Terms generated by expanding the action 5*2 in Eq. ( p. 12 ) to higher than second order in 



ip turn out to be even more irrelevant than ui and U2, as we will see below. We thus have 
a Gaussian FP /i* = (c, 0, 0, . . .) which superficially appears to be stable for d > 4. The 
relevant parameters are t and T or Q, and the correlation length exponent is z/ = l/[t] = 1/2. 
c is marginal, and all other coupling constants are irrelevant. This Gaussian FP obviously 
corresponds to the saddle point solution discussed in the previous subsection. 

The zero-loop RG analysis presented above would imply that the mean-field critical 
behavior is exact for d > 4, unless the RG, at some stage in the renormalization process, 
generated new terms in the action which are relevant for d > 4. From general power counting 
it follows that the only such terms one has to worry about are local terms of 0((5^) in Eq. 
( p.l3| ). Such terms would be relevant and grow like b^. We will see in the next section 



that this indeed happens, because random-field like terms of this type are generated by the 
renormalization process. The net result will be that the upper critical dimension is c?^ = 6 
rather than d'^ = 4, but the mean-field critical behavior is still exact for d > d^ = 6. 
Also, important aspects of the scaling arguments that follow from the above zero-loop RG 
considerations will remain valid. Before we turn to a more sophisticated RG approach, we 
therefore discuss scaling. 

The RG arguments given above imply that the OP satisfies a homogeneity relation, 

N{t,Q,ui,U2,...) = b^"^/^ N{tb^/'',Qb',uib^-'^,U2b^-'^,...) , (3.12a) 

with = 1/2 and z = d/2 + 1. The exponents (3 and z follow from Eq. ( |3.12a|) by means of 
standard arguments.i The crucial point is that u ~ ui ~ U2 is a dangerous irrelevant variable 
(DIV): Solving the saddle point equations explicitly yields N(t,il = 0,m ^ 0) ~ and 
N{t = 0, fi, u ^ 0) ~ u-^/^. Taking this into account, Eq. {^J^ yields /5 = 1/2 and z = 3 
in agreement with Eqs. ( p.9|) . In order to drop the Ui from Eq. (|3.12a|) we have to change 
the scale dimension of N, {d — 2)/2, and the exponent z, to their effective values 1 and z, 
respectively, 

N{t,Q) = b-^ N{tb^/'',Qb') . (3.12b) 

Let us also determine the singular or critical parts Xsing of the thermodynamic suscep- 
tibilities dn/dfj,, 7, and Xs, where 7y = limx-^QCv/T is the specific heat coefficent, and Xs 
is the spin susceptibility. To be specific, let us consider 7y. The singular part of the free 
energy, fsmg, satisfies the relation 

Ung{t, T,u,...) = r ('^+^") /..„,(^6^/^ T¥, ub''-^ ...) . (3.13a) 
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In the critical region, fsing ~ uQ'^ ~ I/m, so that m is a DIV for fsmg as welL This imphes, 

fs^n9{t, T) = r (^+^) fs^ng{tb"''. Tb^ ■ (3.13b) 

Differentiating twice with respect to T yields a homogeneity relation and the critical behavior 

for ^V,sing- 

More generally, we note that all of these thermodynamic susceptibilites scale like an 
inverse volume times a time, so their naive scale dimension is d — z. The DIV changes this 
to 4 — z = 1, so that the result for 'jv^sing actually holds for all of the Xsing, viz., 

Xs^ng{t, T) = b'^^-'^ Xsrng{tb^'\ TV) = ^ ^ Xs^ng{tb\ Tb') , (3.14a) 



or, 



Xsing{t,T = 0) r-. t'/^ , Xs^ng{t = 0,T)r~.T'/'' , (3.14b) 

where we have specialized to the mean-field exponents. 

We next consider transport properties. The charge, spin, or heat diffusion coefficients, 
which we denote collectively by D, all scale like a length squared times a frequency, so that 

D{t, Q) = b'^-' D{tV/\ Vl¥) = b'^ D{tb\ Qb^) , (3.15) 

so the diffusion coefficients scale like the OP, Eq. ( p.l2b| ), in agreement with the explicit 
result, Eqs. (|3.10|) . We are also interested in the scaling behavior of the electrical conduc- 
tivity, a = Dcdn/dfi. The behavior of a depends on whether or not dn/dn has an analytic 
background contribution in addition to the critical contribution given by Eqs. ( p.l4| ). In 
general it is hard to see why dn/d/j,, or any of the thermodynamic susceptibilites, should 
not have an anlytic background contribution, although Ref. |l^ has given some arguments 
why the background may be missing in the particular model under consideration. If dn/ d^x 
has indeed no background term, then 

a(t,^]) = a{tb^,VLb^) . (3.16a) 

The conductivity exponent s, defined by a{t,VL = 0) ~ t*, is 

s = 2u = l . (3.16b) 

If dn/d^ does have an analytic background, then a will scale like the diffusion coefficient, 
so that 



s = 1/ = 1/2 . (3.16c) 

We conclude this section by considering the terms of higher than quadratic order in %l) 
that were neglected in deriving Eq. ( p.l3| ). Expanding the Tr In-term in the action 5*2, Eq. 



( p.lOap , in powers of ip, and localizing each term both in real space and in imaginary time, 
leads to all possible powers of ip. Schematically we denote these terms by 

oo „ 

S^ = Y1 / rfx V^"(x) , (3.17) 

n=l 
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where the Wn are in general matrices. 

Let us examine the Gaussian ip-ip correlation function in momentum space. For c? > 4 it 
goes to a constant at small momenta, for the same reason for which v in Eq. ( p.l5b| ) is finite 



for d > A. If we follow Ref. ^ in defining the scale dimension of a field via the behavior of 



its Gaussian propagator, then we obtain, 

[ip] = d/2 . (3.18) 

We note that other choices for the scale dimension of ip, e.g. [ip] = d — 2 (motivated by 
[Q] = {d — 2)/2 and the fact that A is conjugate to Q^), are also possible, but the final result 
is independent of this choice. 

We next observe that for fixed d, the Wn are divergent if n is large enough. Scaling 
wavenumbers with the correlation length ^, we have Wn ~ ^^n-d ^ y ^^2. Therefore, for 
n > d/2, Wn+i/wn ~ Since each successive term in Eq. ( |3.17] ) carries one more power of 
ip, we obtain for the scale dimensions of the Wn-, 

[w^] = -^(n - 2) for n < d/2 , (3.19a) 



[wn+i] = [wn] - d/2 for n<d/2 , (3.19b) 

and 

d- 4 

[wn+i] = [wn] — for n>d/2 . (3.19c) 

The important conclusions from these power counting arguments are that all of the Wn with 
n > 2 are RG irrelevant for d > 4, that for fixed d they become more irrelevant as n increases, 
and that they all become marginal in c? = 4. Anticipating that we will conclude in the next 
subsection that c/^ = 6, this implies that we were justified in neglecting higher than second 
order terms in ip, provided we work either in d > d'^, or perturbatively (in the sense of an 
e-expansion) in d < d^. 

We conclude this section with three final comments. (1) It is obvious that the nonlocal 
gradient corrections to Eq. (|3.17|) are even more irrelevant than the leading terms studied. 



(2) If one chooses [ip] = d — 2 instead of Eq. ( |3.18| ), then the r.h.s. of Eq. ( |3.19a| ) gets 
replaced by -{n - 2){d - 2) - (rf - 4), and d/2 and {d - 4)/2 in Eqs. ( ^JUBj , |3l9^) get 
replaced hj d — 2 and d — A, respectively. Obviously, all conclusions drawn above remain 
valid. (3) The power counting arguments employed so far do not distinguish between terms 
of different symmetry at a given order in powers of the field. The terms neglected in the 
V'-expansion could lead to terms with different symmetries than the ones present in our Q- 
action, Eq. ( p. 131) , and, because of complications to be discussed in the next subsection, such 
terms could in principle modify the universality classes for the AMT, even though they are 



irrelevant by power counting. We will show in Sec. ^ and in Appendix that all possible 
terms of this kind are generated by the RG anyway, starting with Eq. ( p. 13 ), and that they 
do not modify the critical behavior near d = 6. 
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IV. RENORMALIZATION, AND RANDOM FIELD ASPECTS 



In this section we consider the renormahzation of the field theory derived in Sec. y. First 
we show that a Wilson-type RG procedure generates additional terms in the action. The 
most interesting ones are of random-field (RF) type. Structurally they are identical to the 
terms that would have resulted if the original, unreplicated action had had terms like, 

Srf = ^ j dyi /i„(x) V5^,„(x)^/'^,„(x) , (4.1a) 

with hn a Gaussian random field with a second moment given by, 

{/i„(x)/i™(x)} = e(nm) ^ 5(x - y) . (4.1b) 

Here A/4G' is the strength of the random field, witht the factor 1/46* inserted for conve- 
nience. Standard argumentsEl imply that such an RF term yields c?^ = 6 instead of df^ = 4. 
We then use a 6 — e expansion to describe the AMT in d = 6 — e dimensions. A stable RG 
fixed point is found, and the critical exponents are calculated to 0(e). 



A. Perturbation Theory, and the Renormahzation Group 

Let us consider Eq. (|2.13| ) at criticality, where (A) and (Q) vanish at zero frequency. 
Consider the quartic term, 

^4 = -^ jd^trQ\^) , (4.2) 

and the momentum-shell decompositionj^l 

Q(x) =g<(x)+Q>(x) . (4.3) 

Here has only Fourier components between and unity, with h the RG length scale 
factor, and has only Fourier components between zero and where we have set an 
ultraviolet momentum cutoff equal to unity. Inserting Eq. ( [4. 3D into Eq. ( |4.2| ), examining 
terms ~ (<5^)^, and averaging over the degrees of freedom, yields to 0{u) a term 5*2,2, 

Here ^-^^^ denotes an integral over the region < A; < 1. Examining the structure of 6*2,2 
we see that it is indeed a RF like term, generated by renormahzation at one-loop order. 

Similarly, if we consider terms of 0{v?), we find that renormahzation generates new types 
of quartic terms such as {tr^Q'^){tr_Q'^) and {tr^Q^y. Physically, these new contributions 
refiect the fact that in general all of the terms in the Landau theory of the AMT are random. 
For example, a random quadratic term leads to (tr+Q^)^. All of these quartic terms are RG 
irrelevant for d > 4. However, one of the main technical points about RF problems is that 
the coupling constant. A, of a term like 6*2,2, Eq- ( |4.4|) , scales like 6^ near the Gaussian FP. 
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Since the quartic couphng constants, Ui, scale hke b'^~'^ near this FP, products Aui of a RF 
couphng constant and a quartic couphng constant scale like b^~'^, i.e. the product is marginal 
in d = 6. This leads to d'^ = 6.0 In Appendix ^ we show that, at least to one-loop order, 
of all possible quartic terms only the one that appears in the original action, i.e. utrQ^, 
couples to the RF term. This simplifies the problem substantially. It means that a term 
with the structure of S'2,2, Eq. ( [4.4|) , is the only one that has to be added to Eq. ( p.l3|) in 
order to perform a one-loop RG analysis. Expanding Q about its average value, Eq. ( ^.7] ), 
we can write the effective action as, 



S[4>] 



dx tr 



x)(-9^ + (A))0(x) + |(((g)0(x))^ + (Q)V(x; 



A r 

+77 E {t^ 



u 



2u 



i=>,< 
dx tr 



— Jd^tr [A(0^(x) + -|=(g)0(x 



(Q)03(x 
2 



(ix tr I 



(x)(g)02(x) 



(ix tr 



with A and B given by 

A((Q),(A)) = (g)2-l + /((A)) 



5((Q),(A)) = (A)(g)-2G'/7r] 



(4.5a) 



(4.5b) 



Note that (Q) and (A) are determined by the conditions (0) = {ip) = 0. At zero-loop order, 
these conditions yield A{{Q), (A)) = B{{Q), (A)) = 0. This is just the zero-loop equation of 
state discussed in Sec. pTj 

To perform a loop expansion we need the Gaussian (G) propagator for the field theory 
defined by Eqs. (|4.5|). In the replica limit we find. 



!(k) 

^13(524 + 

with 1 = (ni,ai), etc., and. 



J J,«3«4|'„^l\(G) 



p))(^) = (27r)^ 6{k + p) 



^rs ^ij 



lQ{k'^ + 17112) 



X 



4Ae(nin3) 

C>14^23 H T^— Oro0iodu034 

+ mi2 



(4.6a) 



(4.6b) 



Note that the last term in Eq. ( 4.6a ) diverges like k~'^ at criticality. In propagator language, 
this is what changes the upper critical dimension from c?^ = 4 to c?^ = 6.c3 Alternatively, 
one can leave the RF term out of the propagator and treat it as an interaction vertex.li3 Also 
note that this anomalously divergent critical fluctuation is only present in the spin-singlet 
channel, and only for correlations that are diagonal in frequency and replica space. This is 
in accord with our discussion in Sec. ffl. 



B. The e-Expansion 

We now follow the standard Wilson-type RG approach to construct differential recursion 
relations for the parameters in Eqs. ([4.5|). Let us first consider the critical theory. Then 
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we can put (Q) = (A) = 0, and consider only the renormalizations of A and u. The 
gradient-squared term is not renormahzed at one-loop order, so the field renormalization 
factor, C = is 

( = 5-(d-2)/2+0(62) ^ 

or, equivalently, the critical exponent rj is, 

r] = 0(e2) . (4.8) 
To one-loop order we find the following recursion relations for u and A, 

u' = Ch' \l - lu/\ I ^] , (4.9a) 



2 7k,> A;6 

A' = C^6°' A . (4.9b) 
Taking a differential RG approach, we write 

/ -^^ = ^6 In 6 + 0(e) , (4.10) 

where Sd = Sd/{2TTY with Sd the surface of the {d — l)-sphere. Defining 

g = uA , (4.11a) 

we obtain the fiow equation, 

^^.9 4 s./ + 0(9=) . (4.11b) 

We see that for e > 0, or < 6, there is, besides the unstable Gaussian FP, a new nontrivial 
stable FP. The FP value of g is, 

g* = 2e/9Se + 0{e^) . (4.11c) 



An important feature of Eq. ( [4.9a|) is that the coefficient does not depend on whether 



or not the spin-triplet degrees of freedom {i = 1,2,3) are present. This is different from 
the situation near d = 2, where the presence or absence of the triplet channel changes the 
universality class of the AMT.i Structurally, the insensitivity to the triplet channel in the 
present case is due to the fact that the spin-singlet propagator is more strongly infrared 
divergent than any other term because it couples to the random field. Consequently, the 
leading critical singularities are controlled by the spin-singlet degrees of freedom. 

Next we determine the correlation length exponent, u. Consider the disordered or insu- 
lator phase, where A^„=o = and £„=o = i 0. A straightforward renormalization of (A), 
or i, in Eq. ( [4.5a|) , yields, 

^^=2i-gSd + 0{g') . (4.12) 
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Similarly, a one-loop renormalization of A and B in Eq. ( [4.5a|) gives the equation of state 
to that order, 



N' = l- f{i) 



Gg 



E 



1 



Au V {p^ + 



(4.13a) 



,^ = 2G™-f . (4.13b) 

where both and i should be considered to be functions of fl. Using Eq. (|4.11c| ) in Eq. 
we find, 

In order to determine u, we finally need the relation between £ and the distance from the 
critical point, t. To this end, we expand the r.h.s. of Eq. (|4.13a|) for small i. The i- 
independent contribution is t. At linear order in d = 6 one finds a term proportional to 
i, and one proportional to iln£. The prefactors of these terms are related by Eq. ( |2.15a|) . 
Replacing g hj g*, we can exponentiate and obtain, 

i ^ ti+e/is+o(e^) (4_15) 
Combining Eqs. (4.14) and ( [1.15| ) gives, 

which identifies the exponent i^as l/u = 2 — e/3 + O(e^), or 

1 e . ^ 
z/ = - + — . 4.17 

2 12 ^ ' 

Equations ( [4.8| , [4.17|) give i] and u to first order in e. Standard scaling arguments! yield all 
other static exponents. We find, 

7 = 1 + ^ + 0(6^) , 



/?=l-l + 0(e2) , 
^26 ^ ' 

S = S + e + 0{e^) . (4.18) 

In order to obtain these results we have used the fact that u is dangerously irrelevant even 
for (i < 6, as it is in magnetic RF problems.^ The arguments are the same as those given 
in Sec. for o? > 6, except that Eqs. ( |4.9a| ) and ( [4.11| ) imply that u scales like 

u{b) = u{b = 0) , (4.19) 

instead of u{b) ~ 6^"'^. The net result is that in all ci-dependent scaling laws d is replaced by 
d—2. We note that the exponent values given bvEqs. ([4.8| , |4.17|) and ( [4.18|) are identical with 
the corresponding ones for the RF Ising model.ll3 This may not be too surprising. Since the 
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random potential couples only to the 'longitudinal' field components, oQ°", the problem is 
structurally very similar to an anisotropic RF magnetic model studied by Aharony,@ which 
also yielded RF Ising exponents. 

Finally, we need to determine the dynamical exponent z. For this purpose we consider 
Eq. ( |4.13b| ), which relates A^, and ^2. Expanding the r.h.s. for small £, going to criticality. 



and exponentiating yields. 

Using Eq. ( [4. IS] ) with t replaced by N^l^ yields. 



z = 3-- + 0(6^) . (4.21) 

Note that z = 5f3 /p = yh, with yh the exponent of the field that is conjugate to the OP. This 
was to be expected, since the RG did not couple different frequencies, so that, effectively, 
VL in Eq. ( p. la ) literally acts as the field conjugate to the OP. The technical reason for this 



simplification is that the electron-electron interaction terms turned out to be irrelevant. The 
physical reason is that the static RF fluctuations dominate over the quantum fiuctucations 
which couple different frequencies together. 



V. SCALING DESCRIPTION OF THE ANDERSON-MOTT TRANSITION 

In Sec. |T|we gave scaling considerations that applied to the regime d > c?+, where mean- 
field theory gives the correct critical behavior. Here we generalize this scaling description in 



the light of the RG analysis in Sec. IV 



We first ask which of the general concepts discussed in Sec. ^ will survive. One impor- 
tant general feature is the presence of a dangerous irrelevant variable (DIV), uM Suppose 
that u is characterized by an exponent 9, u{h) ~ . One-loop perturbation theory in Sec. 



IV gave 9 = 2 + 0(e). Although there is reason to believe that ^ = 2 to all orders in the 6 — e 



expansion, this is probably misleading, see the discussion in Sec. and we keep 9 general. 
This adds a third independent exponent to the usual two independent static exponents. In 
addition, there is the dynamical critical exponent z, cf. Eq. (|3.11bD . One effect of the DIV 
is to effectively change z to z. In Sec. |IV| we found z to be not independent, but rather to 
be equal to the exponent yh of the field conjugate to the OP. We saw that this was due to 
the RF fluctuations being stronger than the quantum fluctuations, and the resulting lack 
of frequency mixing. The dominance of RF fluctuations seems to be a general feature of 
RF problems. However, it is possible that in low enough dimensions the explicit interaction 
terms could become relevant, which would lead to frequency mixing, and to an independent 
exponent z. For simplicity, we ignore this possiblity here, and thus write our first scaling 
law as, 

z = y,, = SP/iy . (5.1) 
The OP obeys a scaling or homogeneity relation, 

Ar(t,fi,u,...) = ^(^6^/^^]6^Mr^...) , (5.2a) 
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which upon ehmination of the DIV u turns into 

A^(t, fi) = Ar(t6i/^ . (5.2b) 

This relates the OP exponent j3 to the three independent exponents z/, 77, and 9 through the 
scahng law, 

(3 = ^ i^d-e -2 + 7]) . (5.3a) 

The remaining static exponents are given by the usual scaling laws, with d d — 9 due to 
the violation of hyperscaling by the DIV, 

5 = (^d-d + 2-7])v/2p , 

7 = z/(2 - r^) . (5.3b) 

Now we consider the thermodynamic susceptibilities dn/dn, 7y, and Xs- From the 
general argument given in connection with Eqs. ( p.l4|) we expect all of them to share the same 
critical behavior. Denoting their singular parts collectively by Xsmg again, the generalization 
of Eq. (|37[4il) reads, 

Xsingit, T) = 6-^+^+^ Xs^ng{th''''. T¥) ■ (5.4) 

This links the static critical behavior of the thermodynamic susceptibilities, characterized 
by an exponent k, Xsingif) ~ f^, to that of the OP, 

K = (3 , (5.5) 



where we have used the scaling laws, Eqs. (|5.1| , |5.3b|) . We see that as a consequence of the 



dominant RF fluctuations all of the thermodynamic susceptibilities scale like the OP. This 
is what we had found in mean- field theory, Eq. (|3l4bD , and we now see that this is generally 
valid. 

Finally, we consider the transport coefficients. The homogeneity relation for the diffusion 
coefficients, Eq. ( ^.151 ), does not contain d explicitly, and therefore is generally valid. If we 



denote the static exponent for the diffusion coefficients by sd-, D{t) ~ t"*^, we find the 
following scaling law, 

sd = z — 2 = [3 — PTj , (5.6) 

The mean-field result that the D scale like the OP is therefore valid only to the extent that 
?7 = 0. The behavior of the electrical conductivity cr, which is related to the charge diffusion 
coefficient by means of an Einstein relation, a = Dcdn/dfi, depends again on whether or 
not dn/dfi has an analytic background contribution. Generally one would expect such a 
background to be present. The conductivity then scales like the charge diffusion coefficient. 



a{t,n) = b^-' a{tb^/\T¥) , (5.7a) 
and the static conductivity exponent s (cf. Eqs. (|3.16|) ), is given by 
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s = -(d-2-e-T]) . (5.7b) 

In models or physical situations where dn/dfi has no analytic background one gets instead, 

a{t,Q) = b^-^+^ (x{tb^/'',Q¥) , (5.7c) 

which leads to 

s = v{d-2-e) . (5.7d) 

In either case. Wegner's scaling law s = z/((i — 2) Jll which previously had been found to hold 
for the AMtI as well as the Anderson transitionS, is violated, unless Eq. ( [5.7b| ) holds, and 
6 = 2 — d — f]. This has profound consequences which we will discuss in the next section. 

VI. DISCUSSION 

Our results suggest that a crucial physical aspect has been left out of all previous analyses 
of the AMT in general, and the NLcrM commonly used to describe it in particular. This 
aspect is the random-field (RF) nature of the transition, the presence of which has been made 
plausible in the Introduction, and which has been confirmed by the explicit RG calculation 
in Sec. In our formulation of the problem the RF aspects follow naturally from the 
fact that the random potential couples to the OP for the AMT, cf. Eq. ( [I.2[ ). In this 
sense an OP description of the AMT is necessary in order for the RF features to emerge 
in a straightforward way, and the absence of an OP theory has been the reason why these 
features have not been noted earlier. 

The picture of the AMT that emerges from our OP theory with RF aspects differs 
crucially in many respects, both physical and technical, from the one that had previously 
been obtained by working in the vicinity of d = 2.1 Most importantly, the AMT according 
to the present picture is driven by the behavior of the OP, rather than by the soft modes as 
in the 2 + e expansion. Indeed, we integrated out the soft modes at an early stage in Sec. 
^ and due to the high dimensionality we are working in they never influence the leading 
nonanalytic behavior. Working in a high dimensionality (i.e., in the vicinity of d = 6) was 
in turn forced upon us by the RF aspects of the theory: It is the dominance of the RF 
fluctuations over the quantum fluctuations, which drive the transition in a 2 + e description, 
that shifts the upper critical dimension from rf^ = 4 to (ijt = 6. This is in exact analogy to 
the case of RF magnets, where the RF fluctuations dominate over the thermal fluctuations, 
which is why the RF transition is sometimes refererred to as a zero-temperature fixed point .Ell 
As a consequence of the irrelevance of the soft modes we find only one universality class, 
irrespective of whether or not the spin-triplet channel is present. In contrast, the rich variety 
of universality classes near d = 2 resulted from the influence of the diffusive modes in that 
channel, the number of which can be controlled by adding magnetic impurities, a magnetic 
field, or spin-orbit scattering to the modehili For dimensionalities d < A one expects the soft 
modes to play an important role again. In mean-field theory this is readily seen explicitly, 
see Appendix 0. However, beyond mean-field theory the behavior of the model, and even 
the nature of the transition, are unclear in that regime, as they are in the magnetic models, 
see below. 
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Another important difference between the present picture and all previous work on the 
AMT is that the electron-electron interaction turned out to be irrelevant, in the RG sense, 
for all d > 2. This is another manifestation of the subordinate role that quantum fluctuations 
play in our theory. This leads to an enormous technical simplification over the work near d = 
2, since it eliminates the frequency coupling that made the latter extremely cumbersome.li 
A physical implication of this simplification is that we find orthodox dynamical scaling, 
while near d = 2 there are several critical time scales, and corresponding dynamical critical 
exponents, and one has only what is known as 'weak dynamical scaling'. @i 

Although it is irrelevant for the critical behavior, the electron-electron interaction is 
of course necessary for the AMT to exist, since for noninteracting electrons one has an 
Anderson transition with an uncritical DOS. This point is correctly reflected by the theory, 
as can be seen from Eq. ( p.l4| ) and has been mentioned after Eqs. ( |3.3| ). The electron-electron 
interaction thus plays a rather trivial, although crucial, role in the theory: Though irrelevant 
for the critical behavior, it ensures that the phase transition under consideration is accessible. 
Indeed, for Kg t the critical behavior increases without bound, Gc ^ oo. This suggests 
a number of distinct phase transition scenarios. The simplest one is that for sufficiently 
small interaction constants, or large Gc, the AMT discussed here gets preempted by some 
other phase transition, such as a pure Anderson transition. This scenario is particularly 
likely if = 0, and if the electron-electron interaction is short-rang^ed, since in this case Kg 
is irrelevant near the Anderson transition FP, at least near ci = 2.0 A different possibility 
is that in the above picture the Anderson transition is replaced by an AMT that is related 
to the one studied near d = 2 for the case when either both Kg and Kt are nonzero, or the 
electron-electron interaction is of long range. This picture leaves room for (1) the Anderson 
transition, (2) the type of AMT discussed in this paper, and (3) the type of AMT discussed 
before.li Which transition is actually realized will depend on the relative strengths of the 
disorder and the electron-electron interaction in the various channels. 

All thermodynamic susceptibilities, including dn/dn, are found to be critical in the 
present theory. This is again in contrast to the results in d = 2 + e, where jv and Xs may 
or may not be critical, depending on the universality class, but dn/dn is never critical. The 
reason is that the critical behavior in ci = 2 + e derives from logarithmic singularities in 
d = 2, and dn/dfi, as a frequency integral over a quantity (the DOS) that is only logarith- 
mically singular itself, cannot have any logarithmic corrections to any order in perturbation 
theory.! We note, however, that on general physical grounds one would actually expect a 
nonanalyticity in dn/dn, given a nonanalytic DOS. The frequency integration that takes one 
from the DOS to dn/dn may weaken the nonanalyticity, but it cannot completely remove 
it. In that respect the result in = 2 + e is hard to understand, and the present one is 
physically more plausible. It also bears some resemblance to the case of a Mott transition in 
a Hubbard model, where dn/dfx is also criticahi We emphasize, however, that for the AMT 
a critical dn/dn does not necessarily mean that dn/dn vanishes at the transition, i.e. that 
the system becomes incompressible. As we have mentioned in Sees. |II| and 0, one expects 
in general a nonvanishing analytic background contribution to all thermodynamic quanti- 
ties except for the OP. In special models, or for special parameter values, this background 
contribution could be absent for some or all of the thermodynamic susceptibilities, but one 
would expect a manifest physical reason for this, like, e.g., a symmetry. Ref. |1^ argued that 
the present model might be such a special case. This suggestion was based on an explicit 
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representation of the susceptibilities in terms of the model parameters H, Ks, and KtfB and 
the scaling to zero under RG iterations of the latter, but no more general physical argument 
has been given. Even if the suggestion was correct for this particular model, however, one 
would expect more general models to contain noncritical background contributions. For 
instance, in Ref. ^ it has been argued that local moment effects, which are absent from our 
model, lead to a spin susceptibility that diverges as T — >^ both in the metallic and in the 
insulating phases, but shows no critical behavior at the MIT. Such noncritical contributions 
can enter additively to the critical ones given by Eqs. ( |5.4] ). Of course this would not change 
the critical behavior of anything except possibly the conductivity: Depending on whether 
or not there is a noncritical background contribution to dn/dn^he critical exponent s for 
the conductivity is given by either Eq. (|5.7b|) or by Eq. (|5.7dD .1211 Alternatively, noncritical 
processes not included in our model might lead to bare interaction constants in the action, 
Eq. (p.7|) , that diverge at T = 0. Such terms at most would cause some of the soft modes in 
our model to be absent. Because of the irrelevance of the soft modes mentioned above this 
would not modify the critical behavior for d > 4. 

The remarkable analogy which we have found to exist between the AMT and RF magnetic 
transitions implies that the AMT will also inherit the complications that are known to exist 
for RF magnets, most of which are not quite understood. For instance, it was mentioned 
in Sec. |V| that the exponent 6 is probably equal to 2 to all orders in perturbation theory, 
as it is in RF magnets.S However, it is also known that this result is misleading, and 
that nonperturbative effects are likely to play an important role in general RF problems. 
This is known as the 'dimensional reduction problem' .0 A possible physical consequence 
of nonperturbative effects is the appearance of non-power law dynamical critical behavior, 
and, more generally, features of the RF phase transition problem that resemble those of a 
glass transition.ii If this is true for the magnetic RF problem, then one should expect the 
same for the AMT. A related question is what happens for d < 4. Fisheiil has shown that 
an infinite number of RF-type operators all become marginal in = 4. Together with the 
fact that the soft modes will become important for d < 4 this means that dealing with 
the problem in ci < 4 will require techniques substantially different from the ones employed 
above. 

There are also some problems that are germane to the AMT that have not been ad- 
dressed yet. In this work we have considered only the particle- hole degrees of freedom, the 
particle-particle or Cooper channel has been omitted. In li = 2 -|- e the Cooper channel is 
known to lead to substantial technical problems, which have not been resolved entirely. It 
remains to be seen whether the Cooper channel is amenable to an easier treatment within 
the OP formulation of the AMT. Also, we have restricted ourselves to a short-ranged model 
interaction. Our main motivation for this restriction is that a Coulomb interaction raises 
questions concerning the treatment of screening and the plasmon mode which are most natu- 
rally addressed within the context of the general field theory given by Eqs. ( |1.1| ), rather than 
the NLcrM employed in this paper. Studying the underlying field theory is also attractive 
since the main motivation for using the NLaM, viz. its construction as an effective model 
for the soft modes in the problem, disappears with the realization that the soft modes are 
irrelevant, at least for d > 4. We leave these problems for future investigations. 

Finally, let us discuss experimental imphcations of our results. A result of great impor- 
tance in this respect is the violation of Wegner scaling expressed by Eqs. (|5.7[). This removes 
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the requirement s > 2{d — 2)/d that followed from Wegner scaling (s = v{d — 2V) together 
with the lower bound v > 2/(i.@ This requirement had led to severe problems! with the 
interpretation of a certain class of experiments on doped semiconductors, most notably Si:P, 
which observed s~l/2in(i = 3. While a dangerous irrelevant variable had been envisaged 
for some time as a possible way out of this problem,^ previous attempts to argue for the 
existence of one had been found to be erroneousS 

Another question concerns possible experiments to test whether or not the AMT has 
indeed a RF character. In this context an interesting point is that RF problems are known 
to contain an anomalously divergent correlation function.El The point is that in the presence 
of a random field the quantity {(0(x)) (0(y))}, with O the OP, (. . .) the thermodynamic 
average, and {• • •} the disorder average, is nonzero even if {(0(x))} vanishes. In the present 
case the anomalous correlation takes the form 

C(x,y) = {iV(x)iV(y)}-Ar2 ^ ^g^^^^ 

with A^(x) the unaveraged local DOS at the Fermi level. At the critical point the Fourier 
transform of C behaves like 

C(fc ^ 0) ~ A;-2+'?-^ . (6.1b) 

The same correlation function, but without the anomalously strong divergence due to the RF 
physics, has been considered in connection with mesoscopic systems,^ and it is measurable, 
at least in principle. An experimental observation of the strong divergence predicted by Eq. 
( |6.1b|) would be a strong indication of the presence of RF features at the AMT. It would also 
be interesting to experimentally look for indications of glassy behavior in the form of the 
abovementioned anomalously slow relaxation processes that one expects in a RF system. 
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APPENDIX A: THE MEAN-FIELD SOLUTION FOR D <6 

There is no reason to believe that mean-field theory is a particularly good approximation 
in the physical dimension ci = 3, or anywhere else below d = 6. Here we determine the mean- 
field results for d < 6 anyway. We do so partly for completeness, and partly to illustrate 
certain interesting aspects and problems which at the mean-field level are easily soluble 
because mean-field theory lacks the complications of the RF aspects of the transition, and 
which may be of relevance for a more complete theory in ci = 3. Throughout this appendix 
we suppress the superscript (0) that denoted the mean-field approximation in Sec. |T|. 

Before we turn to the critical behavior we consider the frequency dependence of the OP 
in the metallic phase. For small frequencies we find, 

n\ for d>A , . 

ioT 2 < d < 4 • ^^^^ 



N{n ^ 0) = N{Q = 0) + const x 
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This low-frequency nonanalyticitv, or long-time tail, is the so-called Coulomb anomaly first 
found by Altshuler and Aronovn^ in perturbation theory. Our mean-field theory correctly 
reproduces the perturbative result. In Sec. |rV| we have seen that fl is the field conjugate 
to the OP A^. This allows for an interesting interpretation of the Coulomb anomaly: It 
is the AMT analogue of the fact that in an 0{n) ferromagnet with n > 2 the magnetic 
susceptibility diverges everywhere in the ordered phase.il This is, of course, a consequence 
of the Goldstone modes which become increasingly important for d < A. 

Now we show that the critical behavior is given by Eqs. ( |3.9| ) for all d > 3. We start 
with the equation of state, Eqs. (^.21). It is obvious that Eq. ( |3.5| ) holds independently of the 
dimensionality, so we have (3 = 1/2 for all d. If we expand the equation of state for small 
frequencies at criticality we obtain for d > 3, 

N{Q ^ 0) = . (A2) 

This determines the exponents S = P/uz = 3 for d > 3. For d = 3 there is a logarithmic 
correction to the and for d < 3 the exponent changes. The fact that the exponent is 
(i-independent for d > 3, rather than for d > 4 is a manifestation of the quantum nature of 
the mean-field AMT: A quantum phase transition in d dimensions is expected to be similar 
to the corresponding classical transition in d + 1 dimensions. Note that this is not the case 
for the RF phase transition discussed in the main text because near d = 6 we found that 
the interaction terms, which would lead to frequency mixing, are RG irrelevant. 

In order to determine the remaining exponents we turn to Eq. ( ^.81 ). Obviously rj = for 
all d. In order to determine u we must deal with the divergencies that occur in the coefficients 
u and V in d < A. To protect the divergency, we keep the momentum dependence of u and 
V, which we then scale with the correlation length, ^. The integral for u is dominated by 
the small frequency region, and we find, 

J const for (i > 4 , , 

^ \ const + tV2^4-d for 3 < < 4 " 

The scaling of the correlation length is now obtained by equating the first and the third 
term in Eq. ( p.8|) . We obtain ^ ~ t^^/^, and hence u = 1/2 for all d > 3. Again there are 
logarithmic corrections to scaling in d = 3. Equation ( [A3| ) immediately gives 7 = 2z/ = 1, 
and from z = = 513 / v we obtain z = = 3. The divergence of v is of no consequence: 
Since u is effectively constant for ci > 3, f scales like ^^"'^ ~ ^-{4-^)72^ gQ ^j-^g product vt in 
Eq. ( |3.8|) scales to zero. At the mean-field level, Eqs. (|3.9| ) thus hold for c? > 3. 



APPENDIX B: TERMS GENERATED BY THE RENORMALIZATION GROUP 

In this appendix we show that the new quartic terms generated by the RG, which were 
not in the original action, do not modify the RF fixed point discussed in Sec. at least 
near d = 6 to one-loop order. The argument is that these new terms do not couple to either 
of the coupling constants g or i whose fiow equations are given by Eqs. (|4.11b|) and ( |4.12|) , 
respectively. 

If we repeat the steps below Eq. (|4.2|) to O(u^), then the RF term in the two-point 



propagator, Eq. ( [4.6a|) , leads to the following new quartic terms in the action 
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^4,1 = vi Jd^ tr+ (Q'(x)) tr_ (q2(x 



(Bla) 
(Bib) 



S,_:, = v, Id^tr^iQH^)] [tr+(Q(x))]' , (Blc) 



5*4,4 = J dx [tr+ (Q(x))] " 



(Bid) 



and additional terms given by tr+ in Eqs. ( pT|) replaced by tr_. Note that these quartic terms 
are distinguished from 5*4, Eq. ( [4 .21 ), by the appearance of additional traces. Ultimately these 
additional traces generate extra replica sums that cause these terms not to couple into Eqs. 
( WuBj) and (111). 

First we argue that these terms cannot contribute to Eq. ( [4.12| ). We use Eqs. (^]3|) and 



I]), and examine the (Q^) -terms. The RF singularity is necessary to obtain logarithmi- 



cally singular contributions in = 6, so it is sufficient to use the last term in Eq. ( |4.6a| ) for 
the internal propagator, (Q^Q^)- This either leads to terms that vanish in the replica limit, 
or to terms proportional to [tr±Q]'^. The latter terms are nonsingular renormalizations of 
A, and can be neglected near d = 6. 

Next we consider the renormalizations of 0{uv) and 0(f ^) of S4 or u, Eq. ( |4.2| ), where v 
can stand for any of the coupling constants in Eqs. (pT|). Let us consider the term of ©(fg) 
specifically, the arguments for all other terms follow the same lines. In schematic notation, 
the replica structure of this term is. 



4,2 



Qa\a2 Qa2a\ Qd^ai QaACtd. Ql3i(i2 qP^Pa qPaPz 



(B2) 



We now use Eq. ( [4.3| ), examine the (Q^)''-term, and require that one of the internal prop- 
agators be proportional to the last term in Eq. ( [4.6a|) in order to obtain a logarithmically 
singular term in d = The resulting replica structure is either {tr Q^)^, or {tr QYtr Q^. 
We conclude that this term does not contribute to the renormalization of u. 

Similar considerations lead to the general conclusion that none of the terms given by 
Eqs. (BT) renormalize 6*4, Eq. 
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